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I. INTRODUCTION 

Experimental investigations of the low-energy muon-transfer reactions in collisions of 
muonic hydrogen H^ (bound state of a hydrogen isotope and muon fi^) with nuclei of 
charge Z\ > \ are of importance for muon catalyzed fusion cycle [Q. The study of such 
collisions involving three charged particles is also very interesting from a theoretical point of 
view as an example of rearrangement scattering with Coulomb interaction in the final state. 
Such reactions with post-collision Coulomb interaction between clusters appear frequently in 
atomic and molecular physics [^ . In the following we develop a general formalism for dealing 
with such reactions and as an example apply it to study some muon-transfer processes. 

Recently, there has been considerable experimental interest in the study of the muon- 
transfer reaction in collision of the muonic atoms with He~'"+ |^,^ and also with charges 
2'i > 3 |§-|^, e.g. oxygen (O®), neon (Ne^°), argon (Ar^^) etc. It was found that contrary 
to the smooth Z-dependence expected from the semiclassical Landau-Zener formula [^ the 
experimental muon transfer rates for reactions like 

(H^)i. + X^^X^ + H (1) 

depend in a complicated manner on the charge Z . Here H stands for the hydrogen isotopes 
^H or ^H and X stands for the target nuclei. Another phenomenon which has not yet found 
a satisfactory theoretical explanation is the measured isotope effect, e.g. the trend of the 
direct transition rates of reactions ([1|) for X^ = O^ 0, Ne^° 0, Ar^^ 0, and Xe^^ 0. In 



cases of O*, Ar^^ and Xe^"^ the direct transfer rate decreases with increasing the mass of the H 
isotope. Theoretical analyses |T^ also support this trend. The experimental results for Ne^° 



and sulphur dioxide [§ differ considerably from the theoretical predictions. Moreover, 
several experiments performed in recent years have put into evidence the complex structure 
of the time distributions of the X-rays following transfer from muonic hydrogen isotopes to 



heavier elements 13 



The proper theoretical analysis of charge transfer reaction (|l|) becomes extremely compli- 
cated numerically as the charge Z increases because of the presence of the strong Coulomb 



interaction in the final state. Traditionally, in theoretical studies, such Coulombic systems 
with two heavy (nuclei) and one light (muon) particles are considered within the framework 



of the two-state molecular Born-Oppenheimer approximation |14,15|. In another study, a 



semiclassical model based on Faddeev-type scattering equations has been used [|T6|. It would 
be of interest to perform a full quantum mechanical consideration in view of the fact that the 
muon is not so light compared to the nucleon and compare with the approximate calculations 
mentioned above. 

Here we develop a quantum mechanical approach based on Faddeev-Hahn-type equations 
for a careful reinvestigation of these three-body direct charge-transfer reactions with strong 
Coulomb repulsion in the final state. As a first step towards a model solution of this 
complicated problem, we apply this detailed few-body method to the study of direct muon- 
transfer reaction (|l|) for X^ = ^He''"*', ^He"^'*', ^Li"^"*"^ and ^Li^^^. This study with lighter 
nuclei is expected to lead to faster numerical convergence than the heavier targets. However, 
our approach is equally applicable for heavier targets with higher charges, although the 
convergence could be slow in these cases. These studies with heavier targets would be 
interesting future works. 

For the three-charged-particle system, say (^Li ^H /i), only two asymptotic two-cluster 
configurations are possible, i.e. (^H^) — ^Li and (^Li^) — ^H. For the theoretical treatment 



of such a three-body rearrangement process, Faddeev-type equations [0, especially the 



modified version proposed by Hahn |T8[, appear to be very suitable. The two possible 
asymptotic configurations of the above rearrangement problem are conveniently tackled by 
a set of two coupled Faddeev-Hahn-type equations for components \l/i and \l/2 of the wave 
function \1/ = \I'i -|- \l/2, where each component carrys the asymptotic boundary condition 



for a specific configuration [IS, 20]. These equations are very useful to incorporate distortion 



potentials for specific initial and final asymptotic states [21]. It is possible to include the final 



state Coulomb interaction explicitly in these equations, so that a low-order approximation 



to these equations produces the correct asymptotic behavior ]|2T[ . 

We solve the integro-differential form of the Faddeev-Hahn equation by the close-coupling 
approximation scheme involving up to six states. This procedure consists in expanding the 



wave function components \E'i and ^2 in terms of eigenfunctions of subsystem Hamiltoni- 
ans in initial and final channels, respectively. Although, these subsystem eigenfunctions are 
not orthogonal to each other, the components \l'i and \E'2 satisfy a coupled set of equa- 
tions incorporating the correct asymptotic behavior of the wave function. Consequently, 
there is no problem of overcompleteness as encountered in similar expansion approaches for 
rearrangement reactions based on the Schrodinger equation. The resultant coupled Faddeev- 
Hahn-type equations are then projected on the expansion functions. After a partial-wave 
projection this leads to a set of one-dimensional coupled integro-differential equations for 
the expansion coefficients, which is solved numerically. 

In Sec. II we develop the formalism. We have calculated transfer rates for reaction 
(1) for H = ^H or ^H and X^ = 3He++, '^He++, ^Li+++ or '^11+++ using a two-state close- 
coupling approximation, and for 11 = ^11 and X =^ He''""^, ^Li''""^''' or ^Li"'"'~^ using six-state 
close-coupling approximations. Our results obtained for muon-transfer rates from hydrogen 
to helium and lithium are given in Sec. Ill and compared with those of other investigations. 
We also present a summary and outlook in the concluding part of this section. 

II. THEORETICAL FORMULATION 

Let us take the system of units to he e = h = m^ = 1, where m^ (e) is the muonic mass 
(charge), and denote, the heavy nuclei (^He, '^He, ^Li, etc.) by 1, the hydrogen isotopes 
(^H, ^H or "^H) by 2 and muon by 3. Below the three-body breakup threshold, following 
two-cluster asymptotic configurations are possible in the system 123: (23) — 1 and (13) — 2. 
These two configurations correspond to two distinct physical channels, also denoted by 1 
and 2. These configurations are determined by the Jacobi coordinates (7^-3, /4) 

rj3 = ns - fj, pk = (rg + mjfj)/{l + ruj) - f^, j ^ k = 1, 2, (2) 

fi, rrii are coordinates and masses of the particles i = 1,2, 3, respectively. 

Let us introduce the total three-body wave function as a sum of two components 

^(ri,r2,f3) = ^l(f23,Pl) + ^2(ri3,P2), (3) 



where "^1(^^23, Pi) is quadratically integrable over the variable r23, and "^2(^13, P2) over the 
variable fi^. The components \l/i and \E'2 carry the asymptotic boundary condition for 
channels 1 and 2, respectively. The second component is responsible for pure Coulomb 
interaction in the final state. These components satisfy the following set of two coupled 
equations 

{E-Ho- 1/23)^1(^23, Pi) = (V23 + 1^12 - f/c)*2(ri3, P2) 

(4) 
{E-Ho- Vi3 - f/c)^2(ri3, P2) = (Vis + ^^12)^1(^23, pi) , 

where E is the center-of-mass energy, Hq is the total kinetic energy operator, and Vij{rij) are 
pair-interaction potentials {i ^ j = 1, 2, 3), and Uc is a distortion interaction, e.g. Coulomb 
repulsion in the final state between clusters (^He, //) and ^H in the case of ^He ^H /i system 

Uc^^-^l^^. (5) 

P2 

Here Zi is the charge of ^He and Z2{= 1) is the charge of the hydrogen isotope. By adding 
the two equations (|^) we find that they are equivalent to the Schrodinger equation. For 
energies below the three-body breakup threshold they possess the same advantages as the 
Faddeev equations, since they are formulated for the wave function components with correct 
physical asymptotic behavior. 

The component \E'i carries the asymptotic behavior in elastic and inelastic channels: 



^i(r23,p'i) ^^-^^ e^'^ >i(f23) + Y.^T'^i^py- ''Vn(r23)/Pl . (6) 

The component \E'2 carries the Coulomb asymptotic behavior in the transfer channels: 

^2(ri3,p'2) - E^w(^pJe^^'="''""'^'^''"''^''^'''°''^'''''V™(ri3)/p2, (7) 



ml 



where e 1 ^<^i(r23) is the incident wave, ipn{rjs) the n-th excited bound-state wave function 



of pair (j3), fc» = y2Mi{E - E^J'^), with M"^ = m,^^ + (l + mj)^i . Here E^^^ is the binding 
energy of (j3), i 7^ j = 1,2, A'^/^iVLp^) and A^^iVlp^) are the scattering amplitudes in the 
elastic/inelastic and transfer channels. The Coulomb parameters in the second transfer 
channel are 12! 



n = argr(Z + 1 + V2fc^)) and r] = 2M2{Zi - l)/k^^\ 



(8) 



This approach simphfies the solution procedure and provides the correct asymptotic behavior 
of the solution below the 3-body breakup threshold. 

Let us write down (Hf) in terms of the adopted notations 



E + 



Pk 






^fc(rj3, Pk) = {Vj3 + Vjk - UcSj2)'^j{fk3, Pj) , (9) 



2Mk 2pj 

here j y^ k = 1,2, M^^ = m^^ + (1 + "^j)~^ and pj^ = 1 + 'mj^- We are using the Jacobi 
coordinates 



Pj = ^j3 - Pkns, rj3 = -{Pkpk + Pj) and Vjk = -{(^jpj - crupk) 

7 7 



(10) 



with 



/3fc 



rrik 



cTfc = 1 - /5fc and 7=1- (3k(3j 



(11) 



l + nik 

For solving (y) we expand the wave function components in terms of bound states in initial 
and final channels, and project this equation on these bound states. The expansion of the 
wave function is given by 

1 



LM\l n Pk 



LM ' 



(12) 



where (nlX) = a are quantum numbers of a three-body state and L is the total angular 
momentum of the three-body system obtained by coupling / and A, Vjm's are the spher- 
ical harmonics, Rniy'jz) the radial part of the hydrogen-like bound-state wave function, 
fni\ (Pk) are the unknown expansion coefficients. This prescription is similar to that 
adopted in the close-coupling approximation. After a proper angular momentum projec- 
tion, the set of integro-differential equations for the unknown expansion functions f^\pk) 
can be written as 

92 A(A+1)" 



(1)^2 



dpi 



P\ 



fi'\pi)=9ii: 



(2A + 1)(2A' + 1) 



2L + 1 



Jo 



^ + ^ 



U, 



c 



, , , , K"i\ri3\) 

\r23\ \ri2\ 

P1P2 E DtA0,u,0)C^^^CJ:^;,^,Y;^{u,,n)Y,m'{iy2,r:) 



,(2) 



(13) 



{k^'f + 



dpi 



A(A + 1) 



pi 



-f/, 



c 



fL'\p2)=92T. 



^(2A + 1)(2A'+1) 



2L + 1 



K'l'i\r2s\) 



dpifli\pi) / rf^sm^i?i?(|fi3|) , , , 

insi \ri2 

P2P1 E ^w(0,c^,0)C,^oT™Cf,™:^,y^:„(z/2,vr)l^w(^i, 



Zi Zi 



,{i) 



vr 



(14) 



Here c/^ = AwMk/^f^, 7 = 1 - mfcmj/((l + mfc)(l + m^)), a' = (riTX'), D^^,{0,uj,0) the 
Wigner function, Cxq^^ the Clebsh-Gordon coefficient, cj is the angle between the Jacobi 
coordinates pi and pi', z/j is the angle between rj/3 and pi, u^/ is the angle between r^^ and 
Pit. The following relations are useful for numerical treatment 

Pi' . . 1 



sm Ui 



sin LJ and cos z/,- 



-(APi + Pi'COSu;) (i^i' = l,2). 



(15) 



To find unique solution to (|T^) — ([Tip, appropriate boundary conditions are to be con- 
sidered. First we impose /„; (0) = 0. For the present scattering problem with 1 + (23) as 
the initial state, in the asymptotic region, two solutions to ([T3|) — (p!^) satisfy the following 
boundary conditions 



fil\p^ 



Ji) 



Pi-^+00 



sin(/i;J ' pi) + Kii cos{k[ ' pi 



(1), 



.(2) 



,(2) 



.(2), 



ffKp2) ~ yV^iri2Cos(fcrV2-V2A;r'ln2A;rV2) 



(16) 



where Ky are the appropriate coefficients. For scattering with 2 + (13) as the initial state, 
we have the following conditions 



f(i) 

Is 

c(2) 



fis\pi) ~ Jv2/viK2iCos{k\^'^pi] 



nnP2 



S2), 



12) 



(2). 



P2^+oo 



sin(A;rV2 - v/'^kY' \^2kY' P2) + ^^22 cos(A;rV2- 



J2), 



r^/2kfhYi2kf^P2) , 
where Vi {i = 1, 2) are velocities in channel i. In the absence of Coulomb interaction Uc in 



the final channel, Ky are the components of the on-shell i^-matrix p2[. With the following 



change of variables in ( 13 ) — ( n ) 

ffP(pi) = /L'Hpi)-sin(M'Vi), 

fi?(P2) = fS\p2) - sHkf^p2 - r,/2kf^ \n2kf^p2) , 



(17) 



we obtain two sets of inhomogeneous equations which are solved numerically. The coefficients 
Kij are obtained from the numerical solution of the Faddeev-Hahn-type equations. The cross 
sections are given by 

where i,j = 1,2 refer to the two channels and D = K11K22 — -ft'i2-fGi- When k^ —^0: 
a^^ = (712 ~ 1/^1 • For comparison with experimental low-energy data it is very useful to 
calculate the transfer rates 

^tr = f^tr^^o, (19) 

with V being the relative velocity of the incident fragments and A'^o the liquid-hydrogen 
density chosen here as 4.25 x 10^^ cm~^, because A^j.(A;J — i> 0) ~ const. 

III. NUMERICAL RESULTS 

We employ muonic atomic unit: distances are measured in units of a^, where a^ is 
the radius of muonic hydrogen atom. The integro-differential equations were solved by 
discretizing them into a linear system of equations. The integrals in Eqs. ([T3|) and (|14D 
are discretized using the trapezoidal rule and the partial derivatives are discretized using a 
three-point rule |2^. The discretized equation is subsequently solved by Gauss elimination 



method. As we are concerned with the low-energy limit only the total angular momentum 
L = is taken into account. Even at zero incident energy, the transfer channels are open 
and their wave functions are rapidly oscillating Coulomb waves. In order to get a converged 
solution we needed a large number of discretization points (up to 900) adequately distributed 
between to 40a^. More points are taken near the origin where the interaction potentials are 
large; a smaller number of points are needed at large distances. For example, near the origin 
we took up to 40 equally spaced points per an unit length interval a^, in the intermediate 
region (p = 10 — 20a^) we took up to 25 equally spaced points per unit length interval a^, 
and in the asymptotic region (p = 20 — 40a^) we took up to 15 equally spaced points per 

8 



unit length interval a^. The following mass values are used in the unit of electron mass: 
m(iH) = 1836.152, m(2H) = 3670.481, m(3He) = 5495.882, m{^}ie) = 7294.295, m(*^Li) 
= 10961.892, m(^Li) = 12786.385 and the muon mass is m^ = 206.769. 

We present muon-transfer rates A^^ calculated using the formulation of last section for 
processes (|l]). First, we restrict ourselves to a two-level approximation by choosing in the 
relevant close-coupling expansion the hydrogen-like ground states (H^)is and (X^^Jis, where 
H = ^H and ^H, and X^ = 3He++, ^He++, ^11+++ and ^Li+++. Numerically stable and 
converged results were obtained in these cases. The rates A^^ /lO^ sec~^ at low energies are 



presented in table 1 together with the results of ll^-plq]. The results in this case converged 
to the precision shown in this table, except in the case of ^H^ + ^He'*'''", where it was difficult 
to get converged result. The present results are consistent with the experimentally observed 
isotope effect |l6|,|9|JTT[] , e.g., the rate decreases from ^H to ^H. 



Is 



In table 2 we present our results for transition rate of reaction (|l]) to (^He^^)i5, (^Li 
and (^Li'l^^)is from (^H^)is using the six-state close-coupling model. The six states are 
H^i(ls,2s,2p) and X^(ls,2s,2p). The results so obtained are consistent with the measured 
isotope effect. The effect of including the (2s, 2p) states in the calculational scheme is also 
explicit there. 

The results reported in table 1 and 2 demonstrate the efficiency of the present few- 
body model in describing muon transfer from H isotopes to nuclei of charge Zi = 2. Its 
application to nuclei involving higher charges, therefore, is also expected to be justified. 
The present calculation with ^Li^^^ or ^Li^^^ represents the ffist examples for such a full 
quantum-mechanical extension within the six-state close-coupling model. 

The study of three-body charge transfer reactions with Coulomb repulsion in the final 
state has been the subject of this work. We have studied such reactions employing a detailed 
few-body description of the rearrangement scattering problem by solving the Faddeev-Hahn- 
type equations in coordinate space. To provide correct asymptotic form in the final state 
the pure Coulomb interaction has been incorporated directly into the equations. It is shown 
that within this formalism, the application of a close-coupling-type ansatz leads to satisfac- 
tory results already in low-order approximations for direct muon-transfer reactions between 



hydrogen isotopes and light nuclei He"*"^ and Li"*""*"^. Because of computational difficulties, 
in this preliminary application we have considered up to six states in the expansion scheme 
(ls,2s,2p on each center — (H^) and X^), which may not always be adequate. Further calcu- 
lations with larger basis sets are needed to obtain accurate converged results. However, the 
inclusion of three basis states on each center is expected to build in a satisfactory account of 



the polarization potential in the model. It has been observed [^ in studies of positron and 
positronium scattering using close-coupling type approach that once the ls,2s,2p states of 
positronium and target states are included, a good account of scattering including transfer 
reaction is obtained (estimated error of 10 — 20%). However, the inclusion of only the Is 
basis functions do not lead to the converged results. A similar conclusion can be obtained 



from tables 1 and 2. In view of the results of ref. |^ we do not believe the results of table 
2 to be very different from the converged ones, although we cannot provide a quantitative 
measure of convergence. If the above conclusion based on the works of ref. [^ hold in this 
case we expect a maximum error of 20% in table 2. 

Because of the present promising results for the muon-transfer rates of ([I|) for Z\ < 4, it 
seems useful to make future applications of the present formulation for larger targets with 
Zi > 4. Such calculations involving nuclei of higher charge are in progress. The present 
approach should also be useful in rearrangement collision involving electron, e.g., such as in 
H(ls) + He++ ^ H+ + He+(ls), considered in ^|. 
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TABLES 

Table 1. Low energy muon transfer rates A^^/10^ sec~^ from proton (^H^)is and deuteron 
C'iif,)is to hydrogen-like ground state (^He^)is, ('^He^)is, (''Li^ )i^ and C^Li^ )i^ within 
two-state close-coupling model. 



System 


Energy 

(eV) 


Present Results 

H^(ls),X^(ls) 


1161 


m 


m 


iH^ + 3He++ 


<0.04 
0.1 
1.0 


8.4 
8.3 
8.1 


7.25 


10.9 


6.3 


iH^ + 4He++ 


<0.04 


6.8 


6.65 


10.7 


5.5 


2H^ + 3He++ 


<0.04 
0.1 
1.0 


5.2 
5.1 

4.7 


4.77 


9.6 


1.3 


2H^ + 4He++ 


<0.04 


5.0 ±0.3 


4.17 


9.6 


1.0 


2H^ + 6Li+++ 


<0.04 
0.1 
1.0 


1.2 
1.2 
1.1 


1.01 






2H^ + 7Li+++ 


<0.04 
0.1 
1.0 


1.12 
1.12 
1.06 


0.96 
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Table 2. Low energy muon transfer rates A^j./10® sec ^ from (^H^)is to hydrogen-like ground 

state (^He )is, (^Li )is and C^Li )is within six-state close-coupling model. 

System Energy Present Results 

(e\0 H^(ls,2s,2p),X^(ls,2s,2p) 

2H^-F3He++ < 0.04 9.0±0.2 

0.1 8.8±0.2 

1.0 5.0±0.2 

2H^-h6Li+++ < 0.04 1.9±0.1 

0.1 1.9±0.1 

1.0 1.2±0.1 

2H^-h'^Li+++ < 0.04 1.6±0.1 

0.1 1.6±0.1 

1.0 1.2±0.1 
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